Mirror symmetry breaking as a problem in dynamic critical phenomena 
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The critical properties of the Frank model of spontaneous chiral synthesis are discussed by applying 
results from the field theoretic renormalization group (RG). The long time and long wavelength 
features of this microscopic reaction scheme belong to the same universality class as multi-colored 
directed percolation processes. Thus, the following RG fixed points (FP) govern the critical dynamics 
of the Frank model for d < 4: one unstable FP that corresponds to complete decoupling between the 
two enantiomers, a saddle-point that corresponds to symmetric interspecies coupling, and two stable 
FPs that individually correspond to unidirectional couplings between the two chiral molecules. These 
latter two FPs are associated with the breakdown of mirror or chiral symmetry. In this simplified 
model of molecular synthesis, homochirality is a natural consequence of the intrinsic reaction noise 
in the critical regime, which corresponds to extremely dilute chemical systems. 
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I. INTRODUCTION 

Mirror symmetry is broken in the bioorganic world 
and life as we know it is invariably linked to biological 
homochirality. An outstanding problem associated with 
the origin of life is to explain chiral symmetry breaking 
in nature, why for example, it came to be that the nu- 
cleotide links of RNA and DNA incorporate exclusively 
righthanded sugars while the enzymes involve only the 
lefthanded amino acids. A recent survey of hypotheses 
concerning this phenomenon, experimental realizations 
and additional pertinent bibliography can be found in 
the references [it, 0, i, i, Hi . 

The essential key ingredients of theoretical models of 
mirror-symmetry breaking processes in chemistry |4|| in- 
clude reactions in which the products serve as catalysts 
to produce more of themselves while inhibiting the pro- 
duction of their chiral or mirror-image counterparts. In 
chemistry, enantiomers are molecules that are nonsuper- 
imposable complete mirror images of each other. Frank's 
original model [6^, and a variant of which we study in 
this paper, involves autocatalysis of the two enantiomers, 
denoted here as L and D Q , and mutual inhibition or an- 
tagonistic effects between the two chiral species. More re- 
cently, Sandars introduced a model in which the detailed 
polymerization process and enantiomeric cross-inhibition 
are taken into account, its basic features are explored 
numerically, but without including spatial extent, chiral 
bias or noise Brandenburg and coworkers have ana- 
lyzed further properties of Sandars' model and have pro- 
posed a truncated version including chiral bias 0], and 
have studied this reduction with spatial extent and cou- 
pling to a turbulent advection velocity [l^. Gleiser and 
Thorarinson analyze the reduced Sandars' model with 
spatial extent and coupling to an external white noise 



[ll| and in [l^, Gleiser considers the reduced chiral bi- 
ased model with external noise. Despite the simplicity 
of the Frank model, ignoring as it does the polymeriza- 
tion process, it continues to serve as a kind of "Ising 
model" for chiral symmetry breaking, and the purpose of 
this paper is to better understand its critical properties 
by exploiting the model's relation to directed percolation 
phenomenology. 

The specific reaction scheme we will study in this paper 
is given as follows. The ki denote the reaction rate con- 
stants and we take the achiral substance A as a uniform 
constant background. 
Autocatalytic production: 



L + A^L + L, D + A, D + D. (1) 

h h 

Dimerization and additional mutual inhibition in second 
order reactions: 



L + D 



L + D ^L + A. 



L + D 



D + A. 

(2) 

Spontaneous decay or recycling back to the achiral sub- 
strate: 



L^ A. 



D 



A. 



(3) 



'Electronic address: 
t Electronic address: 



hochbergd@inta.es| 


http : //www. cab . inta.es | 


zorzanomm@inta.es 





The above scheme differs from the original Frank 
model in the open-flow reactor nature of the pro- 
cess and the fact that the reaction Eq.([T]) is allowed 
to be reversible {k^ > 0). We assume that each enan- 
tiomer diffuses with the same diffusion constant Dq and 
incorporate this feature into the master equation descrip- 
tion of this process. We also account for two inhibitory 
or mutually antagonistic reactions, with associated rates 
ki and k^, in addition to Frank's dimerization step, fc2. 
This scheme is a partial hybrid between the Frank model 
and the Avetisov and Goldanskii (AG) reaction (see e.g., 
Eq.(13) of [H). Whereas Frank's original model gives rise 
to pure homoquiral states in which only one enantiomer 
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is present, the complete AG model leads to chiral sym- 
metric broken final states where both enantiomers are 
present in unequal proportions. Mirror symmetry is bro- 
ken in the AG model, but the breaking is not absolute. 

Both the above and Frank's 1953 scheme yield the 
same field-theoretic structure for the effective action, and 
more importantly, therefore belong to the same univer- 
sality class. In our analysis, we allow for ^ fcs, as 
this leads to a rich fixed point structure in the critical 
regime of the model. This choice implicitly accounts for 
the infiuence of an external chiral field or bias. Of course, 
a chiral symmetric action and Langevin equations result 
from the "natural choice" k4 = k^. The properties of the 
chirally unbiased model can be studied as a special case 
of the above scheme. 

Chiral or mirror symmetry breaking is an example of a 
nonequilibrium phase transition which is attained when 
the control variable {kiA — kg) « becomes small. A 
continuous, second order, transition is then induced from 
a fully active state, characterized by the simultaneous 
presence of both competing enantiomers accompanied by 
fiuctuations in each chirality, to an inactive or absorbing 
state, in which only one enantiomer survives. This con- 
trol parameter can be made small (or large) by simply ad- 
justing the concentration of the achiral molecule A. This 
limit implies that the net amount of total chiral matter is 
vanishingly small, i.e., the chemical system is extremely 
dilute at criticality. This is because the autocatalytic 
amplification of the enantiomers is delicately balanced 
by their spontaneous decay. The purpose of this paper 
is to understand this specific critical behavior and the 
emergent properties of the Frank model as exposed quan- 
titatively by applying results [13] from the field-theoretic 
renormalization group (RG). We describe a significant 
result in the field of molecular chirality, namely, sym- 
metry breaking induced by internal reaction noise in ex- 
tremely diluted systems with absolute enantioselective 
catalysis. Such a mechanism is vitally important for cur- 
rent scenarios of prebiotic chemistry, where it is com- 
monly accepted that sufficiently high concentrations of 
organic compounds could not have been reached during 
the chemical evolution of the early earth. 

We are interested in the long time and long wavelength 
properties as governed by the nature of the RG fixed 
points and the topology of the RG flow in the space 
of effective reaction rates. The statistical field theory 
derived from the scheme Ea. (|ll2l3p maps identically to 
an action for so-called multi-species directed percolation 
(MDP) Ts], for the special case of two "colors" or species. 
Thanks to this correspondence, the full details of the RG 
analysis already carried out for MDP can be carried over 
and applied to analyze the critical chemical properties 
of the Frank model. In the next Section, we present 
the field-theory action associated with the above scheme, 
which after a suitable rescaling, leads to the effective ac- 
tion that holds in the critical regime. In Section IIIII we 
reproduce the complete RG flow diagram for this model. 
However, only a part of this flow diagram is applicable to 



real chemical systems, and we discuss the consequences 
for chiral symmetry breaking near criticality. In Section 
HVl we derive the Langevin equations that individually 
hold in the vicinity of the saddle point and the two stable 
fixed points of that fiow diagram and integrate these nu- 
merically to obtain the time dependence of the competing 
enantiomers for both large and small noise amplitudes. 
The results are briefly summarized in SectionlVlwhere the 
signiflcance of criticality for scenarios of prebiotic chem- 
istry is emphasized. The relation between criticality and 
extremely diluted chemical systems is brought out in Ap- 
pendix [X] The modiflcations that must be made to the 
effective action when the dimerization and antagonistic 
reactions are allowed to be reversible are briefly discussed 
in Appendix [Bl 



II. THE EFFECTIVE FIELD THEORY ACTION 

The correct inclusion of the effects of microscopic den- 
sity fluctuations in reaction-diffusion systems can be car- 
ried out once the kinetic scheme is specified. With the 
scheme in hand, we derive the corresponding chemical 
master equation, represent this process by creation and 
annihilation operators on a spatial lattice [14], and in the 
final step, upon taking the continuum limit, we pass to a 
path integral representation [l^[l^. From this, an effec- 
tive action Seff can be straightforwardly derived which 
contains all the critical dynamics implied by the reaction 
scheme to be studied. The mapping of related kinetic 
schemes to continuum statistical path integrals is spelled 
out in [13, [3 where the main steps can be found. Ap- 
plying this procedure to the scheme in Eqs.lH)),® and 
([3]) yields the complete action S governing the reaction 
dynamics: 

S = J d'^xj dt |a* (^dta - DaV^a + kiab - kiAa 
+ fcgfl + k^a^^ + a*'^ (kj,a^ — kiAa) 
+ b* (dtb - DoV^b + fcaafo - kiAb + k^b + k^b^^ 

+ b*^{k3b^ - kiAb) + k2a*b*ab + kib*ab + kia*b*ab 
+ k5a*ab + k5a*b*ab} , (4) 

where d is the spatial dimension, and 
a(x, i), a*(x, i), 6(x, i) and 6*(x, t) are continuous 
fields. In the absence of noise (the mean field approx- 
imation) the fields a(x, t), 6(x, t) correspond to the 
coarse-grained local densities of the L and D enan- 
tiomers, respectively. With the noise properly restored, 
these fields are generally complex-as is the noise-and do 
not directly represent the physical densities. However, 
the spatial averages (a(x, i)), (6(x, t)) are indeed real 
and do correspond to the particle densities [l9|. The 
quantities a*(x, t), &*(x, i) represent the conjugate or 
response fields. These are intimately related to the 
fluctuations inherent in the system. In fact, when 
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the action S depends quadratically on the conjugate 
fields, these can be integrated out exactly from the 
path integral, and the noise statistics completely and 
rigorously characterized p^ . The non-critical spatial 
dynamics (i.e., for {kiA — ke) >> 0) implied by the 
action Eq. fl]) with its attendant complex noise and fields 
was explored numerically in for A:4 = /cs = and 
A:6 = 0. 

The fields are next rescaled in the action Eq.(|4]), for 
the purpose of determining which couplings (i.e., which 
combinations of the rate constants ki) are going to be ir- 
relevant in the strict sense of the RG. This step is needed 
in order to correctly identify the complete set of vertices 
that are required to construct a field-theoretic perturba- 
tion expansion of this action . We make use of the 
observation that when k2 = and k^ = k^ = 0, the ac- 
tion Eq.Q reduces to that for two identical uncoupled 
copies of the single particle Gribov process. We there- 
fore rescale the fields according to a* = 9ip*, a — 0~^tp, 
b* = Ocj}* and b = Q-^cj), where 9 = (^)^/^. The 
space and time dependent densities of L ancl D are given 
by ijj and 0, respectively. We define the new coupling 
uq = (fci^fcs)^/^. Introducing a length scale and 
measuring time in units of (i.e., [Dq] = we find 
that the new fields have scaling dimension k''/^, while 
[r] = [{ke — fciA)/_Do] = is a relevant perturbation 
in the RG sense. On the other hand, [uq] = so 
this nonlinearity becomes marginal in dc = 4 dimensions. 
Thus, we learn that d = 4 is the upper critical dimension 
of the Frank model, below which the mean field approxi- 
mation is incorrect. Note that [k2] = [fca] = [^4] = [k^] = 
K^^'', and hence these couplings are irrelevant compared 
to [uo]: indeed, since e.g., [^2/^0] = k~'^^^, these par- 
ticular rate constants may be omitted from the effective 
action. Doing so, leads to the effective action, which 
takes the form: 



S. 



d'^x / dt {^*[dt+Do{r-V^)]ilj 



- wo(V*V-V'>') 

+ 4>*[dt+Do{r~W^)]cl>~uo{cf>*^^-cl>*cl>'') 

+ (fc2 + fc4)0~ + ik2 + k^)e~^i^*u] ■ 



(5) 



Note that the decoupling of the two enantiomers occurs 
when ^2 = fc4 = fcs = 0, in other words, for vanishing 
mutual inhibition and dimerization. In particular, we 
see that both the original Frank model and the extension 
treated in this paper do indeed lead to the same field- 
dependent structure for the effective action, Eq®. 



III. 



CRITICAL BEHAVIOR 



The field theoretic renormalization group (RG) can be 
applied to Sef / in order to study the nonequilibrium crit- 
ical dynamics of this reaction-diffusion system (for a ped- 
agogical review of this methodology, see [l6t ) . The main 



purpose for employing RG techniques is that they lead to 
differential equations describing how the model parame- 
ters, in this case, the kinetic constants, transform under a 
change of length scale. As we are here interested in the in- 
frared, or long wavelength, properties of the Frank model, 
we therefore consider the RG flow of the parameters in 
the long wavelength limit. In general, certain combina- 
tions of the kinetic constants will flow to various flxed 
point (FP) values that depend on the space dimension. 
Thus, the flow diagram can be constructed revealing the 
critical properties of the underlying kinetic scheme. As 
it turns out, Seff maps exactly to a field theory of so- 
called multi-species directed percolation (MDP), for the 
special case of two colors or species [l^l • A complete and 
exhaustive RG analysis has already been carried out for 
the general model in [l3|, and as pointed out there, the 
required renormalization factors for MDP are provided 
by the single species Gribov or directed percolation pro- 
cess. As an immediate consequence, the parameter com- 
bination [kiAk^y-^"^ / Dq = uo/Do flows under renormal- 
ization to the stable fixed point u* — i-y/2e/3, where 
e = 4 - d > 0. 

We next turn to the two interspecies couplings, which 
from Eq. ([5]) are each seen to be proportional to the sum of 
the rates U12 c>c (/c2+^5) and U21 oc (fc2+^4), respectively. 
The competition between the two enantiomers comes in 
through the dependence on the rate of dimerization k2, 
as well as through fc4 and ^5. The complete RG analysis 
in [13] as applied to our model proves that, except for the 
point D, the interspecies parameters in Seff will flow to 
one of the following d-dependent fixed point values: 



(ui2,U2l) 



({k2+k5) 



(fc2 + ki 



Do 
D: (0,0), 
Ul : (0,2m*; 



Do 
S : {u\u*) 
U2 : (2u*,0) 



(6) 



The point D corresponds to complete decoupling between 
the two enantiomers, S to a chiral symmetric coupling, 
whereas Ul and U2 each correspond to homochiral final 
states. The flow of the interspecies couplings U12 and 1*21 
under renormalization is depicted in the flow diagram in 
Figure [1] The flow, as indicated there by the sense of 
the arrows, corresponds to the critical large wavelength 
and long time properties of the microscopic model defined 
in Eqs.([T]), ^ and ([3]), and is reached for small values 
of r « 0. Thus the system goes critical when the differ- 
ence in the rates of autocatalytic amplification {kiA) and 
spontaneous decay (fcg) goes to zero. This is achieved 
by varying the concentration of the achiral matter A. 
This corresponds exactly to a situation of extremely di- 
lute net chiral material characterized by + ~ (see 
Appendix [X] for a simple proof of this fact). Note the 
topology of the fiow and the stability property of each of 
the fixed points (D,S,U1,U2): (totally unstable, saddle 
point, stable, and stable), respectively. 

We first consider the fiow properties as depicted within 
the positive shaded quadrant U12 > and U21 > in 
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FIG. 1: (Color online) RG flow or "phase diagram" of the 
reaction scheme. Space of scale dependent renormalized cou- 
plings Ui2 and 1121. The shaded quadrant corresponds to 
the critical properties of the model with positive kinetic con- 
stants and diffusion. There is an unstable flxed point D at 
(mi2,«2i) = (0,0) (blue dot); a saddle point S at (Mi2,ii2i) ~ 
(m*,u*) (red dot), and two stable fixed points Ul and U2, 
located at (1*12, W2i) = (0,2it*) and (ui2,W2i) = (2«*,0), re- 
spectively (green dots). The separatrix 1121 = 1*12 is the 
boundary between the basins of attraction of Ul and U2. 
Outside this quadrant, the line M21 -|- U12 = separates the 
complete basin of attraction of Ul and U2 from the dashed 
line of fixed points U12 + U21 = —2u*. See text and reference 
[1^ for further details. 



Fig. [TJ In so far as it is reasonable to assume that the 
kinetic constants ki and the diffusion Dq are nonnegative 
parameters, this the most pertineiit part of the full flow 
diagram for real chemical systems [2(j| . There is a totally 
unstable fixed point D located at the origin = "21 ~ 
0. This corresponds to complete decoupling of the two 
enantiomers. But the only way to arrive at this decoupled 
state is by blithely setting the initial values of A:2 = as 
well as fc4 — — all to zero. Otherwise, the slightest 
positive deviation of any one of these rates from zero, 
drives the system eventually to either the saddle point 
S, if dui2 — Su2i > 0, or to one of the two stable fixed 
points Ul, if Su2i > Sui2 > or U2, otherwise. 

The diagonal symmetry line U21 = U12 in this quadrant 
is a separatrix dividing the basins of attraction of the two 
unidirectional fixed points. For chirally symmetric kinet- 
ics, the natural choice of course is k^ — k^, which puts 
the system dynamics directly on top of this separatrix. 
Then, as the diagram indicates, any positive initial value 
for U21 = U12 > drives the system to the chiral sym- 
metric fixed point S. In this case, the final state of the 
system is determined by the fully symmetric couplings 



between the two enantiomers (see Eq.(l6l)). In chemistry, 
a racemic mixture is one that contains equal amounts of 
left- and right-handed enantiomers of a chiral molecule. 
When the system is near the point S, racemic initial con- 
ditions lead to a racemic final state, while non racemic 
initial conditions lead to a final state that maintains the 
original enantiomeric excess only for low noise amplitudes 
(see the numerical results in Section lIVp . 

On the other hand, if the the model has /c4 — fcs ^ 
which from Eq.(l6]) implies that U21 7^ M12, then the sys- 
tem evolves to one of the two stable fixed points Ul or 
U2. At either Ul or U2, the system attains unidirec- 
tional interspecies couplings, see Eq.®, which lead to 
the absolute amplification of one enantiomer at the ex- 
pense of the other: that is, complete chiral symmetry 
breaking and a pure homochiral stable final state is the 
inevitable outcome (see the numerical results in Section 

EYJ. 

For the sake of completeness, we now address the re- 
mainder of the flow diagram (the unshaded regions). In 
this case, there is then another separatrix whose equa- 
tion is M12 + U21 = 0, which divides the complete basin 
of attraction of the two unidirectional fixed points from 
the dashed hue uu + U21 = —2m*; see Figure [T] The 
parameter domain to the left of this dashed line corre- 
sponds to a region of instability, and it is conjectured in 
jJl] that couplings satisfying the condition uu + U21 < 
will lead to first order transitions. We hasten to point 
out, however, that this part of the diagram is only acces- 
sible if initial values of the U12 and or U21 are negative, 
corresponding to a negative diffusion Dq < 0, provided, 
of course, that none of the reaction rates ki are allowed 
to become negative '20]. Thus, the region of this diagram 
applicable to real chemical systems is represented by the 
shaded quadrant. 



IV. CRITICAL DYNAMICS 

The temporal evolution of the two enantiomers in the 
critical regime represented in FiglT] is governed by a pair 
of coupled Langevin equations which follow straightfor- 
wardly from the effective action Seff- These are obtained 
by carrying out a Gaussian integration over the conjugate 
fields ■0* and (f>* in the path integral of the exponentiated 
effective action: J VipV'tp*V(t)V(j}* e'^""^'^''^' ''^''^'l This 
final step yields a product of delta functional constraints 
under the integral which in turn, lead to a pair of ex- 
act coupled stochastic partial differential equations [l^. 
The advantage of obtaining the Langevin equations in 
this way is that the noise properties are fully determined 
and do not have to be guessed at or put in by hand. 
Numerical solutions of these stochastic equations can be 
carried out to reveal the nature and qualitative tendency 
of the spatial and temporal evolution of the competing 
enantiomers in the neighborhood of each RG fixed point, 
as well as within their respective basins of attraction. 

The Langevin equations that follow from Seff are 
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given by 

d 

~ {k2 + h)0-^U + ii (7) 

d 

— = 1)0 + [kiA - ke)(t) - UQcf)^ 

- (fc2 + A;4)e~W^ + 6, (8) 
where the noise satisfies (^i) = (^2) =0 and 

(ei(x,i)a(x',i')> = 2«oV'(x,0'5''(x-x')<5(t-t'),(9) 
(e2(x,i)6(x',t')> = 2Mo0(x,t)(5''(x-x')<5(t-i')(lO) 

These equations hold in the critical region shown in Fig. 

m 

A. Langevin equations in the vicinity of the saddle 

point 

In particular, the behavior of the model in the vicinity 
of the saddle point S where the interspecies couplings U12 
and U21 flow to a symmetric fixed value, is given by the 
solutions of the system Eqs. (|7l8p where we now set ^4 = 
k^. These are subject to the noise Eas. (|9ll0p . and we 
use the result that uq/Dq flows to u*, together with the 
corresponding fixed point values for U12 and 7i2i, as given 
in Eq.®. At this juncture, it is also convenient to rescale 
the fields -ip = Dou*/{kiA - ke)^, = Dou*/{kiA - 
kQ)(j), and employ dimensionless time r = {kiA—kQ)t and 
coordinates xj = {{kiA — k^)/ D^Y/'^Xj. These simple 
steps yield the stochastic equations in the vicinity of the 
saddle point S: 

= V'^^ + ij-^'^ -^^ + Ii{^,t), (II) 

= + -H + U^,t), (12) 

where the rescaled noise is given by 

(er(x,.)^;(x',r')) = 

X 5'^(x-x')'5(r-T'), (13) 

(6(X,r)e2(X',r')) = 2(^^A^) " '^(*'^) 

X J''(x-x')(5(t-t'). (14) 

In two dimensions, the noise strength is characterized 
by the parameter ct^ = 2Dou*^ / {kiA — fee), with u* = 
I/V3 w 0.58, and this can be large or small depending on 
whether the diffusion rate Dq is large or small (keeping 
the difference kiA — ke > fixed), respectively. In Fig- 
ure m some representative effects of the diffusion on the 
critical dynamics in the neighborhood of the saddle point 
S in d = 2 are displayed for the spatially averaged enan- 
tiomers for both large (cr = I) and small (cr = 0.3) inter- 
nal noises and for non-racemic initial conditions. Recall 



that here we have set ki = k^. We solve numerically the 
full stochastic two-dimensional version of Eqs. pilI2p 
subject to the noise given by Eqs. (|I3II4p , using reflect- 
ing boundary conditions and a finite difference scheme 
with At = 0.005, Ax — Ay — 0.23, and a grid of size 
L X L = 154 X 154. 




I ^ ^ 1 

10000 20000 30000 

T 

FIG. 2: (Color online) = k^: Time evolution of the 
spatially averaged enantiomer densities (u) = {(p) and of 
(v) — (■(/)) in two dimensions d = 2 for two representative sim- 
ulations of the stochastic dynamics near the saddle point S, 
Eas. (|llll2p with noise Eqs. (|13ll4p . Evolution of non-racemic 
homogeneous initial conditions (^{cj>{x,y,T = Q),il){x,y,T = 
0)) = (0.6,0.4). For weak noise ct = 0.3 see the inner pair of 
red (dark gray) and green (light gray) curves. For stronger 
a = 1 noise, see the outer pair of blue (dark gray) and ma- 
genta (light gray) curves. 

For small noise levels (cr = 0.3) , and for non-racemic 
initial compositions, the initial proportion of the two chi- 
ral species is roughly maintained, modulo the fluctua- 
tions; see the inner pair of curves in Fig. [51 The evo- 
lution of (u)cr=o.3 is shown in the red (dark gray) line 
and (w)o-=o.3 in the green (light gray) line. However, 
for stronger noise (cr = I), the initial imbalance shows 
an almost monotonic tendency to increase, suggesting 
that sufficiently strong noise is capable of driving the sys- 
tem to a homochiral final state, in spite of the manifest 
mathematical chiral symmetry of the underlying evolu- 
tion equations and noise terms under the substitutions 
tp <j) and (j) ^ ip. A mean field analysis of the solutions 
of Eqs. pTjl and (|I2p. which ignores both diffusion and 
noise, indicates that the enantiomeric excess of the con- 
centrations of the two enantiomers is time independent 
[2l| . Here, this is seen to be approximately true also for 
the spatially averaged diffusing enantiomers subject to 
small noise. But greater internal noise induces a striking 
departure from this that is not captured by the mean field 
approximation, as seen in Fig. [2j This is depicted in the 
outer pair of curves. The evolution of {u)a-=i is shown 
in the blue (dark gray) line and {v)^^i in the magenta 
(light gray) line. 
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As we are interested here in displaying only the initial 
and intermediate time dependent tendencies of the two 
enantiomer densities in the vicinities of the various RG 
fixed points, we have employed standard integration of 
the Langevin equations, sufficient for revealing the qual- 
itative nature of the solutions for short and intermediate 
computational time steps, as can be seen in Fig. 2 and 
Fig. 3 (below). Near an absorbing state transition, one 
of the densities tends to zero, and the numerical inte- 
gration breaks down. This can be seen clearly in the 
simulation of the evolution to the absorbing state, which 
is shown only for the shorter time scales in these figures. 
This standard algorithm is of course not adequate for ex- 
tracting the much more precise and delicate information 
such as asymptotic decays or power law exponents. For 
the latter, we would have had to appeal to the more so- 
phisticated numerical schemes such as those pro posed by 
Dickman [23|, Moro [H or by Dornic et al. [2l|. 



B. 



Langevin equations in the vicinity of the 
unidirectional fixed points 



The behavior of the system near one of the two stable 
attracting fixed points, for instance Ul, is determined by 
the pair of equations 



-2V'</) + a(x,T), (15) 
-6(x,t), (16) 



with the noise properties as given above in Eqs. (|13ll4p . 
Here, we use the fixed point value (mi2,U2i) = (0, 2w*). 
Recall in order to arrive at this fixed point, we set 
^4 7^ fcs . This corresponds to "starting" the system ofi^ in 
either the basin of attraction of Ul or that of U2 (see 
shaded quadrant in Fig. [1]). Note the manifest asymme- 
try in the equation pair due to the presence of the unidi- 
rectional coupling term in Eq. (|15p . absent from Eg. pop . 
This fixed point is associated with homochirality, as con- 
firmed by numerical simulation; see Fig. [31 Starting 
from racemic initial conditions, the plot of the spatially 
averaged enantiomeric densities, in Fig. [31 indicates an 
extremely rapid onset of absolute chiral amplification. 

The evolution of {u)a-=i is shown in the upper solid 
red (dark gray) curve and that of {v)cr=i is shown in the 
lower dashed green (light gray) curve. The equations 
governing the critical dynamics at the other stable fixed 
point U2 are had by simply interchanging the fields V' ^ 
(j) in Eqs (|15ll6[) and replacing i^i ^ ^2- The behavior 
is quantitatively identical, but with the roles of the two 
enantiomers obviously reversed. 

Near criticality, there is no numerical evidence for 
the formation of the spatially segregrated chiral domains 
bounded by racemic fronts, in marked contrast to the re- 
sults reported in [l^ . The distinguishing factor of course 
is that the present simulations are carried out at the uni- 
directional critical points Ul or U2, whereas in [l^, the 
system was explored far away from criticality, where such 




FIG. 3: (Color online) ^ kz,: Time evolution of the 
spatially averaged enantiomer densities {u) = {<j)) and of 
{v) = (-i/j) in two dimensions d — 2 for a representative simula- 
tion of the stochastic dynamics near the asymmetric unidirec- 
tional fixed point Ul, Eqs. p5ll6|l . Evolution of homogeneous 
racemic initial conditions ((0(£,y,r — 0),ijj(x,y,r = 0)) = 
(0.5,0.5) for noise level a — 1. {u)a=i is shown in the upper 
solid red (dark gray) curve and {v)a=i is shown in the lower 
dashed green (light gray) curve. 



chirally pure domains are expected to form [y, [2^. At 
the RG unidirectional fixed points, the dynamical equa- 
tions themselves are manifestly chirally asymmetric, and 
the system rapidly evolves to a final homochiral state, 
without passing through the intermediate stages of enan- 
tiomeric competition. 



Stochastic dynamics away from the critical 
points 



The distinction between the critical and non-critical 
behavior can be sharpened by contrasting the mathe- 
matical structure of the critical equations and noise to 
those that hold away from the fixed points. The latter 
are of course given by the system Eas. (|7l8p with the fluc- 
tuations obeying Eas. (|9ll0p [26]. By means of the field 
rescaling ^ ~ {{k2 + k^)/uo)'ilj, and with a similar rela- 
tion between and ((>, these Langevin equations can be 
written as (note: in what follows we take = k^) 

dr-tjj = V^^ + -0-g^/i^ -^/i0-|-?7i(x, r), (17) 

dr4> = \/''4> + 4>~94^^-4^^ + V2{i,T), (18) 

where g = kz/(k2 + k^) and the rescaled noise obeys 
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(r?i(x,r)77i(x',r')) = 2'^!^^^{k,A ~ hY/^-^{1.^T) 

X 5''(x-x')<5(t-t'), (19) 
(,72(X,r)r;2(x',r')> = 2^-^^^±^{hA - hf/^-'^'i^^r) 

X 5'^(x-x')'5(T-r'). (20) 

If the noise and the diffusion terms are ignored, then a 
mean field analysis of the homogeneous asymptotic solu- 
tions of the corresponding equations Eas. (ll7ll8l) . reveals 
that the parameter g plays a special role [2Il- Indeed, 
f; < 1 leads to chiral amplification of an initial enan- 
tiomeric excess, whereas t; > 1 leads to a racemic fi- 
nal state. The point g = 1 was identified as a critical 
value, in the sense that a sudden qualitative change in 
the asymptotic behavior of the mean- field solutions is ob- 
served: the ratio of the two enantiomeric concentrations 
in this borderline case remains constant, and depends on 
the initial composition. If the initial condition is racemic, 
the system will always remain racemic, if however there 
is a slight initial excess, this excess is forever maintained 

From these remarks we see that the spatially de- 
pendent and stochastic effective equations associated 
with each renormalization group fixed point S,U1, 
Eas. (|llll2p and Eqs. pSllGp . respectively, have g — I- 
Likewise for the point U2. It is as if we had set g to its 
"critical" value. But it is important to emphasize that at 
these RG fixed points, g is no longer a freely adjustable 
parameter, but under renormalization is automatically 
driven to this special value. The RG thus provides a 
rational physical explanation for why g = 1 at criticality. 

V. DISCUSSION 

As supported by surveys and reviews of the present 
status of chiral autocatalysis, mirror symmetry breaking, 
stochasticity, and their implications for the origin of ho- 
mochirality, the Frank model and its extensions continue 
to serve as the paradigm for theoretical studies of this 
phenomenon [j, H, [23| . In this paper we have studied 
the critical properties of the Frank model and a simple 
extension of it, by exploiting the mapping of this kinetic 
scheme to a well studied phenomenon from condensed 
matter and non-equilibrium statistical physics, namely, 
(multi-species) directed percolation processes [l^. By 
virtue of this exact mapping, which is established at the 
level of statistical field theory, the complete renormaliza- 
tion group (RG) analysis of the critical properties of di- 
rected percolation can be applied to study the critical fea- 
tures of the Frank model. The most significant result in 
this paper is the description of a new effect in the field of 
molecular chirality, namely mirror symmetry breaking in- 
duced by internal noise in extremely diluted systems with 



absolute enantioselective catalysis. Such a mechanism 
is of course vitally important for scenarios of prebiotic 
chemistry, where it is commonly agreed upon that suffi- 
ciently high concentrations of organic compounds could 
not have been reached during the chemical evolution of 
the early earth [28(]. These final states are consequence 
of internal composition fluctuations and reactions limited 
by spatial diffusion. To reach these dilute multi-critical 
states, the difference between the amplification and de- 
cay rates must be close to zero [l^. This contrasts to 
Saito and Hyuga's suggestion that, for closed systems, 
both nonlinear autocatalysis and recycling with diffu- 
sion seem to be req uired for chiral symmetry breaking in 
dilute solutions [2^. Additional insight into the dynam- 
ical consequence of each transition is provided by deriv- 
ing and numerically solving the exact effective Langevin 
equations that hold in the neighborhood of each fixed 
point. 

It has been known for some time that the chirally sym- 
metric state of the Frank model is unstable and that (ex- 
ternal) fluctuations can induce a transition to homochi- 
ral final states. There are evidently a number of distinct 
routes leading to homochirality, and the concept of criti- 
cality and the identification of the associated critical pa- 
rameters should be clearly distinguished. From general 
bifurcation theory, we thus learn that the transition from 
a symmetric to a chiral final state can be induced by vary- 
ing solely the concentrations of the substrate molecules 
[3g |. The mean field analysis of Ref. [2l| on the other 
hand, identifies the ratio of rate constants g = k^/k2 as 
the pertinent critical parameter. In certain crystalliza- 
tion experiments, it is the stirring rate of the solution 
that has been observed to play the role of a critical pa- 
rameter [3l|. The present work makes use of the fact 
that chiral or mirror symmetry breaking is an example 
of an active to absorbing state phase transition, and that 
such transitions are generically characterized by directed 
percolation processes (DP)[i^. RG techniques can be 
applied to analyze this symmetry breaking phenomenon 
in extremely dilute chemical systems in a precise manner. 
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APPENDIX A: CRITICAL REGIME IMPLIES 
EXTREME DILUTION 

Simple bifurcation analysis of the purely kinetic 
scheme Eas. (|ll2l3p 

^ = {kiA~ke)L~k3L'' -{k2 + k)LD, (Al) 

^ = ikiA-ke)D-k3D^ -{k2 + k)LD, (A2) 

reveals the chemical nature of the critical regime of the 
fully stochastic field theory treated in this paper [s^ ]. 
Here, A, L and D denote concentrations and we have set 
k4 = k^ = k. Introduce the enantiomeric excess rj = 
and the total concentration of chiral matter x ^ L + D. 
Then the kinetic equations Eqs. (|AllA2p can be written 
as follows: 



drj k2 + k — ks 



(A3) 



dt 2 

— = {kiA-k6)x-[k3 + ^ (1-'7)]X- 

(A4) 

For behavior with large amount of chiral matter, x » Oi 
we must have {kiA — fcg) >> 0. This requires that the 
system be far from criticality. The bifurcation equation 
[l| is then 77(1 — 77^) = and there are three stationary 
solutions: 



{V = 0, X = 2 }, 
k2 + k + ks 

unstable if fc2 — fc > and stable if A;2 — A; < 0, and 



(A5) 



{V = ±1,X = ; 1, 



(A6) 



stable if fc2 — > and unstable if fc2 — fc < 0. From this 
we can deduce the following salient features: First, by 
introducing the dimensionless time r = {kiA — kQ)t, we 
see that far from criticality, the variable x changes much 
more rapidly than the enantiomeric excess 77. The system 
rapidly reaches a quasistationary state for x {dx/dr w 0) 
and then the slow variable 77 evolves and the full system 
reaches its true steady state. Secondly, the criticality 
condition {kiA — fee) ~ corresponds to the kinetic be- 
havior under conditions of extreme dilution, x ~ 0, with 
the concentration of chiral material close to zero. In this 
case, the system has no well defined steady state with 
respect to the concentration x, since 77 = [L — D)/x 
yields an indeterminate expression near criticality. Fur- 
thermore, the equation for x becomes as "slow" as the 
equation for 77, and the use of classical kinetic approach 
based on the law of mass action becomes questionable. 



Thus, the stochastic approach employed in this paper is 
not only justified, but is needed to correctly describe the 
critical regime of the Frank model. 

APPENDIX B: REVERSIBLE REACTIONS 

The backreaction of the dimerization step is eliminated 
if the product P is continuously being removed from the 
reactor. Otherwise, the reverse reaction must be taken 
into account. Furthermore, as pointed out by Avetisov 
and Goldanskii reversibility in the mutual inhibition 
reactions will account for the limited cnantioselectivity of 
chiral catalysts, so that the catalytic effect of each enan- 
tiomer leads to the formation of both L and D products. 
To include these effects, the following reactions would 
have to be added to the above scheme Eas. pi2l3p : 



k-2 



L + A 



L + D, 



D + A 



fe-r, 



L + D. 
(Bl) 

The lefthandmost corresponds to the backreaction of the 
dimerization, while the latter two allow for reversibility in 
the mutual inhibition reactions. Going through the same 
algebraic procedure that led us to the effective action in 
Secini we find that the above reactions yield the following 
terms to be added to the effective action in Eq.([5]): 



AS', 



d". 



dt{k^2d->P* + k-2i 



+ fc_2eV*r + fc-40*1^ + fc-5?/'* 

+ k^A0ip*^l}(t)*k_z6(j3*(j)i}*} . 



(B2) 



From dimensional analysis we find that [fc_26'] = k^^''^^. 



[fc_20'] = [fc-4] = [fc-5] 



are relevant perturba- 



tions in the sense of the RG for all dimensions, whereas 
[fc_46'] = [fc-5^] = K^"''/^. The corresponding cubic 
terms are therefore marginal in c?c = 4 dimensions and 
are relevant for d < 4. 

From the point of view of the field-theoretic content 
of ASeff, we see that the dimerization backreaction in- 
duces new relevant terms not present in the original ef- 
fective action proportional to '--^ ^*,~ 0*, as well as a 
term of dimension k^. However, this additional term 
proportional to t/)*^* dynamically couples the two enan- 
tiomers via a cross-correlated noise, a feature not present 
in the absence of dimer breakup. Regarding the lim- 
ited cnantioselectivity reactions, these induce "masslike" 
terms ~ that also serve to link the two 

enantiomers. In fact, these terms lead to off diagonal 
contributions to the response functions or propagators. 
The new cubic terms ip*ip(j)* , (j)*(f)'ijj* are also sources 
of "off-diagonal" or cross-correlated reaction noise; in 
graphical perturbation theory, these lead to new cubic 
vertices which would have to be included in a field- 



theoretic RG analysis, such as in [17 1 
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